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ABSTRACT 



A coupled oceanic-atmospheric boundary layer model which 
provides single-staticn prediction capability is evaluated 
relative to wind stress and oceanic mixed layer depth obser- 
vations. The model is initialized and verified using data 
obtained during the 1983 Mixed Layer Dynamics Experiment 
(MILDEX ) . Model computation of friction velocity (u*) and 
oceanic mixed layer depth are compared with observations 
from both atmospheric frontal and non-frontal synoptic situ- 
ations. Favorable results of model u* predictions are 
achieved although in some cases they are slightly higher 
than observed turbulence derived values of u*. Mixed layer 
predictions are very close to the observed except in frontal 
regions. Proper AB1/0BL model initialization of boundary 
layer values is critical in order to achieve favorable 
results. 
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I. INTRODUCTION 



Naval oceanographers and meterologists obtain a great 
variety of information to help them in making atmospheric 
and oceanic forecasts. Extensive use of large shore based 
computers, remote sensing equipment and world-wide communi- 
cations networks make the forecaster rely on remotely based 
technology rather than local indicators and have, to a great 
degree, removed the "art" from single-station forecasting. 

However, large-scale conflict or operational considera- 
tions could deny all forecasting inputs available to the 
forecaster except local indications. Reliance on single- 
station forecasting for naval operations during times of 
crisis or war was first addressed by Oliver and Oliver 
( 1945 ) : 

During the last few years, wartime conditions have made 
it necessary for isolated combat units to issue fore- 
casts in regions where no network of meteorological 
stations could be available. Frequently the data from 
several stations or from reconnaissance planes are 
available, but in seme regions the forecaster must rely 
only on surface and upper-air observations made at his 
own station. This is particularly true in the case of 
ships at sea. Hence, it is important to develop profi- 
ciency in extracting information from limited aerolo- 
gical data. 

While the above quote pertains to atmospheric forecasts, the 
same could be true for ocean predictions. 

In modern warfare, isolated ships or battle groups far 
from the main force, operating in an emission control 
(EMCON) posture, must be able to predict their operating 
environment without shore support. Such predictions are 
important for the effective utilization of weapon and sensor 
systems and to deny the effective use of the environment to 
the opposing force. Affected systems include; navigation 



radars, communications equipment, over-the-horizon (OTH) 
radars, air and surface search radars, missile and gun fire 
control radars, laser and infra-red (IE) guided missiles, 
hull mounted sonars, towed arrays, and variable depth sonars 
(VDS) . 
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SOUND VELOCITY 
STRUCTURE 

ACOUSTIC PROPA- 
GATION MODELS 



Fig. 1.1 Components of the Coupled Atmospheric 
and Oceanic Boundary Layer Model. 



A coupled atmospheric and oceanic model (O’Loughlin, 
1982) has been developed to predict ocean and atmospheric 
parameters from local indications. The components of this 
coupled model are depicted in Fig. 1.1. This model is being 
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tested and evaluated with shipboard applications in mind. 
Accurate predictions of ocean mixed layer depths ( MLD) , 
thermocline jumps and gradients are essential in optimizing 
anti-submarine warfare (ASH) systems and employing acoustic 
countermeasures. These parameters are dependent on the heat 
fluxes and wind stress at the ocean surface. This thesis 
will examine how the single-station coupled model 
predictions of the surface fluxes and mixed layer depth 
compare with observations. 
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II. MODE! BACKGROUND ANE ANALYSIS PROCED URES 



A. BOUNDARY LAYER DESCRIPTION 

Thermodynamically and dynamically forced marine atmos- 
pheric (MABL) and oceanic (OBI) boundary layers must be 
considered in coupled model description. The air-sea inter- 
face has adjacent oceanic and atmospheric turbulent mixed 
layers which effectively insulate the guasi-geostrophic 
oceanic and atmospheric flows above the top of the MABL and 
below the bottom of the OBL. The primary source of the 

turbulence is the velocity (shear) and buoyancy (density) 

gradients created by the exchange of mass, energy and 

momentum across the air-sea interface. Such turbulence 
generates turbulent kinetic energy (TKE) and vigorous 

vertical mixing (Davidson and Garwood, 1984). 

Under undisturbed conditions, the MABL consists of a 
moist, well-mixed layer extending from the surface to the 
capping inversion. At the inversion, jumps in temperature 
and humidity occur with a rapid increase in temperature and 
a corresponding decrease in humidity. Within the MABL, 

equivalent potential temperature and specific humidity are 
conservative quantities. 

A consequence of the large velocity fluctuations and 
mixing in the MABL is that energetic eddies (extending from 
the surface to the inversion) entrain warm, dry air and 
bring momentum into the mixed layer. This results in the 
upward growth of the layer. Subsidence arises from large- 
scale atmospheric forcing and tends to limit the depth of 
the mixed layer. Therefore, the change in the inversion 
height with respect to time is a function of the entrainment 
rate and subsidence. If entrainment causes the MABL to 
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extend above the lifting condensation level (LCL) , clouds or 
fog will form within the layer and greatly affect the OBL. 

In the ocean, the OBL is defined as extending from the 
surface to the top of the seasonal thermocline. Within the 
mixed layer, nearly homogeneous profiles of temperature and 
density are evident. While density is a function of both 
salinity and temperature, over short time scales, tempera- 
ture has a far greater effect (Biller, 1976). At the bottom 
of the mixed layer, large changes in temperature and 
salinity exist with increasing depth. This is the top of 
the thermocline region and serves as a transition from the 
mixed layer to the usually dynamically stable interior 
ocean. Idealized model profiles for both the KABI and OBL 
are depicted in Fig. 2.1. 

Radiation is another source of energy in both mixed 
layers. While most solar energy will penetrate the MA3L, 
this is not true of the OBL. Long wave radiation is 
absorbed and emitted by the first few millimeters of sea 
water. Therefore downward turbulent heat flux is as impor- 
tant as the upward flux. Surface heat flux in the OBL is a 
function of sensible and latent heat as well as the long and 
short wave solar radiation. This radiation is highly depen- 
dent on cloud formation in the ABL so proper ABL modelling 
is essential for realistic ocean MLD predictions. 

While buoyancy forces create turbulence in both boundary 
layers, relatively more turbulence in the OBL is produced by 
wind stress (Davidson and Garwood, 1984) . The MLD will 
deepen if the wind generated downward flux of momentum and 
forced convective mixing (of surface water that has been 
cooled by radiative heat loss or evaporation) provide suffi- 
cient TKZ to erode the thermocline and entrain cooler water 
downward. If, on the other hard, the near surface absorp- 
tion of solar radiation produces a downward buoyancy flux 
which dampens turbulence and provides for net warming, the 
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Fig. 2.1 Idealized Atmospheric and Oceanic 
Boundary layer Temperature Profile. 

MLD will shallow. Fig. 2.2 shews the importance of mechan- 
ical mixing to the CBL and its interaction with buoyancy 
fluxes. A primary objective of this thesis will be to 
determine the effects of wind stress on the MLD. 

Several cause and effect relationships are evident. 
First, clouds can be caused by changes in the ocean surface 
temperatures which, in turn, affect the radiation budget of 
the OBL which results in surface temperature changes. 
Secondly, the stability of the MABL is influenced by air- 
ocean temperature differences. Again, cloud formation plays 
an important role. Such relationships make a coupled 
prediction approach necessary. This coupled model uses 
local oceanic descriptions from a bulk oceanic mixed layer 
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WORK 8Y 
WiNO STRESS 



U.b. 

NET SURFACE BUOYANCY FLUX 
PLUS ABSORPTION Of SOLAR RaOIATKX 




Fig. 2.2 Hechanical Energy Budget for the 
Ocean flixed Layer. 

model and the atmospheric descriptions come from an atmos- 
pheric boundary layer model. Such a coupled model has 
tactically significant implications to forecast parameters 
that are important to weapons systems. Fig. 2.3 shows the 
interrelationships between the coupled model and tactical 
products which can use model outputs (on the extreme right 
hand side) . 
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MABIBE ATMOSPHERIC BOUNDARY LATER (MABL) MODEL 



Existing integrated MABL nodels are zero-order, two 
layer, mixed layer models consisting of a well mixed turbu- 
lent boundary layer underneath a relatively non- turbulent 
free atmosphere. A transitior zone (inversion) of zero 
thickness (hence the term zerc-order) separates the two 
layers. Because it has no thickness, a jump occurs at the 
inversion of the conservative parameters. The model was 
described by Stage and Businger (1981), who formulated 
entrainment energetics, and modified by Davidson, et al. 
(1984) , who described bulk aerodynamic formulae for the 
surface layer. 

The following inputs are required by the model: 

(1) an initial atmospheric sounding; 

(2) the mean winds at a level within the surface layer; 

(3) the surface temperature; and 

(4) the subsidence rate. 

As formulated, ten winds can be input over the forecast 
period. The sea-surface temperature (SST) remains unchanged 
and the surface current is assumed to be zero. SST and the 

surface current change in the coupled version due to the 
interaction with the OBL model. Inversion height, mixed 
layer values of temperature, humidity and wind are predicted 
at 30-minute intervals in both the coupled and uncoupled 
versions. The prediction steps are shown in Fig. 2.4. 
Procedures are the same for clear and cloudy cases except 
for entrainment computation and estimation of cloud top 
cooling. Formation of clouds cr fog, cloud top cooling and 
associated entrainment are important in the physical 
processes in the MABL. 
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Fig. 2.4 Schematic of Input, Prescription and 
Computing Steps in E&BL Prediction. 

The MABL model is very sensitive to the prescribed large 
scale subsidence which moderates the depth of the moist 
marine layer and hence, cloud formation. Methods which can 
be used to compute the subsidence from single station obser- 
vations are: the kinematic method; the adiabatic method; 
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and the integration of the moisture budget equation (Q) 
method. Gleason (1982) found that the Q-method has the most 
merit and was used in this thesis. However, this method 
assumes that no advection will take place. 

Integrated rate equations by Tennekes and Dreidonks 
(1981) are used to predict the time rate of change of the 
conservative quantities and their respective jumps at the 
inversion. The vertically integrated equations for a 
conserved property, X, are: 



h (DX/Dt) = (w’x’) - (w'xM + source (2.1) 

o h 

h (DAX/Dt) = hrx(3h/3t) - (w*T») + (v»xM - source (2.2) 

o h 

where : 

source = - (F -F }/pC D for X = temperature 

source = 0 for X = humidity 

and Gamma ( r ) is the conservative property lapse rate above 
the inversion and F n is the net radiative flux. The 
subscripts "h n and "o" refer to inversion height and surface 
values respectively. 

The entrainment velocity parameterization is based on 
the formulation by Stage and Businger (198 1) and is used to 
close the system of equations and determine the time evolu- 
tion of the inversion height. The closure assumpton is that 
the dissipation rate of the TKE is a fixed fraction ( 1 — A) of 
the production rate. "A" is the entrainment coefficient and 
is taken to be 0.2. 

Bulk aerodynamic formulas described by Davidson, et al. 
(1984), are used to determine the surface fluxes of 
momentum, sensible heat and latent heat. They are: 

u* = c d ^ 2l I 10 (friction velocity) (2.3) 



2 1 



T* = 2(9 0 - 


9) 


(temperature parameter) 


(2.4) 


= <y/2(g 0 - 


S) 


(moisture parameter) 


(2.5) 


These fluxes are 


given 


by: 




u * w * = u* 2 




(momentum) 


(2.6) 


T « v • = 




(sensible heat) 


(2.7) 


g • v * = u*g* 




(latent heat) 


(2.8) 



where C and C. are t*he stability dependent wind and temper- 
ature drag coef f icie nts, 9 is the potential temperature and 
g is the specific humidity. The subscript "o" denotes the 
surface value. 

Because radiation is a primary factor in the OBI devel- 
opment, considerable effort has gone into calculating the 
radiation budget in the MAB1. Uncertainty in background 
aerosols, distribution and concentration of atmospheric 
absorbing gases, and cloud droplet size spectra are sources 
of error in the MABL radiation calcula tions. The short and 
long wave fluxes are computed separately. Short wave radia- 
tive flux is calculated using the delta-Edding ton method 
(Joseph, et al., 1976) which accounts for heating of the 
mixed layer by solar radiation. Incident flux at the top of 
the mixed layer is obtained from the flux at the top of the 
atmosphere and the average transmittance in each of 15 bands 
covering the spectrum from 0 . 2 urn to 1.7um. Short wave 
extinction is due to both scattering and absorption. Water 
vapor in the atmosphere is the primary absorbing constit- 
uent. Scattered radiation forms a second short wave radia- 
tive component (diffuse solar radiation) which, in the mixed 
layer, is due to atmospheric particles (cloud droplets and 
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aerosols). Total short wave radiative flux at the surface 
is determined by computing both the direct and diffuse radi- 
ative components. Fairall, et al. (1981) gives a comprehen- 
sive review of this method. Solar zenith angle is necessary 
in order to compute short wave radiative flux and it is 
determined from the latitude, julian day and time of day of 
the initial input. The fraction of short wave radiation 
from the sea surface, Ag, which is an important short wave 
ratiative parameter, is prescriled in the MABL model as 0.1. 

long wave radiative flux is calculated differently for 
the cloudy and cloud free cases. In the cloudy case, the 
only clouds permitted are non-black stratus. Radiative flux 
is a function of cloud emissivity, e c , which depends on the 
total liquid water content, W, cf the cloud. W profiles are 
nearly linear with height above the LCL (Davidson, et al. , 
1984) and, along with e c can be described by: 



where p g is the density of air (1.25 X 10“ 3 gm/cm 3 ), h is 
the height of the mixed layer (cloud top), Z c is the ICL 
(cloud bottom), a = 0.158 (Slinco, et al., 1982), and q,. is 
the liquid water content at the cloud top. 

Using the Stef an-Eoltzmann law, the net long wave cloud 
top radiation flux, L nh , can be calculated from the cloud 
top temperature, T h , and the effective radiative sky 
temperature, T sky , cloud bottom net radiative flux, L nc , is 
calculated in the same way using the cloud bottom tempera- 
ture, T c , and the SST, T g . These fluxes are given by: 




(2.9) 



€ c = 1 - exp (-aR) 



( 2 . 10 ) 




( 2 . 11 ) 




( 2 . 12 ) 



23 



where a is Stefan’s constant (4.61 X 1 0 _ 1 1 ) and e c is 

obtained from equation 2.10 . Ihe net long wave flux at the 
surface, F, , becomes: 

long ' 



r 

long 



= or [T £ 



4 - 



€ T 4 

c 



- (i - 



6 c) ^Sky 4 ] 



(2. 13) 



where T is the average temperature of the cloud. 

For the cloud free case, the net fluxes are calculated 
at the top of the mixed layer (Z = h) and at the surface (Z 
= 0) by integrating the flux emissivity profile (Fleagle and 
Businger, 1980) , which is a function of the water vapor and 
temperature profiles. The net long wave flux at the surface 
for the clear sky case is: 



long 



= F u- 



(2. 14) 



where F u and F d are the upward and downward radiative fluxes 
respectively. 



C. OCEANIC BOtJNDABY LAY EH (OBI) MODEL 

The OBI model is a one-dimensional, second order turbu- 
lence closure, vertically integrated (bulk) model of the 
upper ocean surface turbulent boundary layer or mixed layer 
developed by Garwood (1977) . It uses the continuity equa- 
tion for an incompressible fluid, the first law of thermody- 
namics (heat equation), the conservation of salt equation, 
an analytical equation of state, the Na vier-Stokes equation 
of motion with the geostrophic component eliminated, and a 
two-component vertically integrated TKE budget. 

An understanding of the dynamics of the entrainment 
process is necessary to predict the MID change with time. 
The stable water mass underneath the mixed layer is destabi- 
lized and eroded by TKE in the nixed layer. This TKE budget 
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is the basis for an entrainment hypothesis. A closed system 
of equations is obtained by csing the bulk buoyancy and 
momentum equations with the mean turbulent field modeling of 
the vertically integrated equations for the individual TKE 
components. 

To better define the mixing process, separate vertical 
and horizontal equations for TFE are used. 3uoyancy flux 
and shear production provide energy for vertical mixing, 
with buoyancy flux being a somewhat more efficient source of 
energy for mixing due to its direct contribution to the 
vertical component of TKE. The buoyancy equation is derived 
from the heat and salt equations using a linearized equation 
of state: 



P = P G [ 1 - « (9 - 9 0 ) - 0 (S - S c ) ] (2. 15) 

Buoyancy is given by: 

b = g (p Q - p)/p Q (2. 16) 

where : 

9 = temperature 
S = salinity 
g = gravity 
p = density 

« = thermal expansion coefficient 

/3 = density coefficient for salt 

The tilde represents the total instantaneous value and the 
subscript "o” denotes an arbitrary, but representative, 
constant value. In the short term, salinity has little 
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effect on the upper ocean density profile except at higher 
latitiudes (Miller/ 1976) . This leaves temperature as the 
dominant density factor. However/ by using buoyancy (b) 
instead of only temperature (9) allows the model tc be 
applied in cases where evaporation and/or precipitation 
contribute significantly to the surface buoyancy flux. 

Model initialization requires: 

(1) mixed layer temperature and salinity profiles; 

(2) wind-driven horizontal current profiles. 

The initial salinity profile is made isohaline if unknown 
and the initial currents are assumed to be zero if unknown, 
lack of initial information about current and salinity 
profiles is not a serious deficiency because the model will 
evolve reasonable transient profiles of salinity and 
momentum after only a single diurnal cycle/ and the final 
results are not sensitive tc the initial current and 
salinity profiles (Davidson and Garwood, 1 984). It is also 
possible to prescribe an upper ocean internal vertical 
velocity (upwelling and downwelling) if it is known to be 
significant. The mixed layer depth, h, is defined as the 
shallowest depth at which the observed density value, cr t , is 
0.02 a t units greater than the observed surface density 
value. 

At each one hour model time step, the following boundary 
conditions are required: 

• net upward turbulent heat flux at the water surface 
(sensible and latent heat flux) plus back radiation; 

• incident solar radiation ; 

• the fraction of short wave radiation absorbed in the top 
one meter of the ocean; 

• surface wind speed and direction; 
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• clcud cover ; 



• temperature (SST, dry bult and dew-point air tempera- 
ture) ; and 

• precipitation (P). 

Vertical velocity at the bottom of the mixed layer may be 
prescribed hourly. Other physical and model constants that 
need to be prescribed are: 

• the extinction coefficient for short wave radiation 
absorption ; 

• the Coriolis parameter; 

• the critical Richardson number for stability adjustment 
below the mixed layer; 

• expansion coefficient for temperature; and 

• density coefficient for salt. 

The model forecast is not particularly sensitive to these 
constants, but reasonable values have been determined by 
Gallacher, et al. (1983), and they need not be readjusted 
for geographic and seasonal variability. 

By use of bulk aerodynamic formulas, the turbulent 
fluxes of latent heat, Q e , and sensible heat, Q h , can be 
estimated as follows: 

°e = C d <-98 E s - V O 10 (2.17) 

0 h = C .( T , - ( 2 - 18 > 

Net back radiation is estimated from the empirical equation 
(Husby and Seckel, 1 578): 
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(2. 19) 



Q„ = 1.14 X 10-7(273. 16 + T ) 4 ( . 3 9 - .5E Vz) 

b S3 

(1 - . 6 C 2 ) 

where: 

E s = saturated vapor pressure of the marine air (0.98 
corrects for salt defects) 

E a = vapor pressure of air based on dew point temperature 

T g = air temperature 

T s = sea surface temperature 

C = fractional cloud cover 

C d = drag coefficient 

The upward heat flux, Q u , is then given as; 

Q u = Q e + Q h + Q to (2. 20) 

and the solar radiation, Q s , by: 

Q s = (1 - a a b ) (1 - . 66C 3 ) Q 0 (2.21) 



The constants "a" and "b" are adapted from Tabata (1964) and 
the cubic cloud cover correction from Laevastu (1960). 
Clear sky radiation, Q 0 , is given by Seckel and Beaudry 
(1973) : 



Q 0 = A q + A, cos $ + B-, sin 0 + A 2 cos 20 
+ 3 2 sin 20 



( 2 . 22 ) 



where the coefficients (A 0 , A 1 , etc.) were calculated by 
harmonic representation of the values predicted in the 
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Smithsonian Meteorological Tables (List 1958) with: 

(J) = (2V365) (t - 2 1) (2. 23) 

where " t" is the julian day of the year. 

The ocean is not very transparent to solar radiation. 
In the open ocean, approximately 50% is absorbed within the 
first meter (this fraction is called "Rf"). Absorption 
varies from region to region and is highly dependent on the 
concentration of absorbing particles such as phytoplankton, 
"Gelbstoff " (yellow substance) , and suspended particulate 
matter. Coastal regions absorb more radiation because of 
increased amounts of suspended particulates. Very little of 
this absorbed radiation penetrates below the mixed layer 
because cf the low thermal conductivity of the underlying 
stable thermocline. Thus most cf this energy is transferred 
upward out of the ocean and back into the atmosphere. What 
short wave radiation remains penetrates the mixed layer and 
is attenuated in an exponential fashion depending on water 
turbidity. Net heat flux at the surface is given as: 

Q net * £ u + ( R f)Qs - Qs (2.24) 

Surface buoyancy (heat and salt) and momentum fluxes can 
be computed by use of the foregoing equations. Mixed layer 
turbulent temperature, salinity, velocity and buoyancy 
fluxes are given by: 



(T ’ w ' ) 


^net /P ^ p 


(2.25) 


(S' w’) 


= (P - E) S Q 


(2. 26) 


(u ’ w') 


= u* 2 


(2.27) 


(b ' w ' ) 


= g[« (T'w») - /3 (S' w') ] 


(2. 28) 



29 



where the subscript "o" refers to the surface value. The 
friction velocity in air, u*, is calculated by: 

r s = Pa C d U 10 2 < 2 - 29 ) 

u* = ( fg /p a ) 0 (2. 30) 

where r g is the surface stress (dynes/cm 2 ) . The fluxes of 
momentum, radiation, latent and sensible heat at the sea 
surface determine shallowing (retreat) or deepening of the 
mixed layer by entrainment. If there is a positive buoyancy 
flux (Q t <0 and E>?) the MLD will deepen. This will be the 
case at night or whenever long wave radiative cooling plus 
the upward turbulent fluxes of heat and moisture exceed 
solar radiation. Negative buoyancy flux will result in the 
shallowing of the MID and is caused by the domination of 
daytime solar heating (net warming) at the surface, provided 
wind stress is not large. 

A schematic of the ocean model computations is shown in 
Pig. 2.5. New ocean temperature, salinity and wind-driven 
current profiles are predicted at one-hour intervals and are 
used to predict the new MLD. 

D. COUPLED BOUNDAEI LAYER MODEI 

Coupling of the atmospheric and oceanic models was first 
accomplished by O’Loughlin (1S82) by matching momentum, 
sensible heat, latent heat, and radiation at the air-sea 
interface. The atmospheric part of the coupled model 
computes the radiative, heat acd momentum fluxes which are 
inputs to the ocean model. Since the ocean model has an 
hour time step, and the atmospheric model has a half-hour 
time step, the ocean model was inserted as a subroutine into 
the atmospheric model and is called every other atmospheric 
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Fig. 2.5 Schematic of Input, Prescription and 
Computing Steps in OBL Prediction. 

time step. Before the atmospheric fluxes can be utilized by 
the ocean model, they must be converted to the appropriate 
units for seawater. 

The ocean model prediction requires the vector surface 
stress for the ocean turbulent velocity flux, u* w 2 . This 
requires that the wind field be broken into horizontal 
components: 



(2.31) 



°10x = ~^ sin 9)U 10 
°10y = " (COS 9) U 10 {2 ‘ 32) 

where U 1Q is the wind speed at 10 meters an ^ 0lO* • ^10y are 

the east-west and north-south horizontal components respec- 
tively. The wind direction, 0, is relative to true north. 
Then u* w 2 is given by: 

T s = p a u* 2 = P w u* w 2 (2.33) 

the density of seawater, p w , to be about 1 gm/cm 3 . Finally, 
the horizontal components of the turbulent velocity fluxes, 
u* wx 2 and u* wy 2 , are computed by: 



u 



*2 = 



Vio 



2 = C,(U 



Q 1 10x 



2 + U 



lOy 



(2. 34) 



u* 2 = 2 

wx Q lOx 



(2. 35) 



U * 2 — 

u wy 






(2. 36) 



SST is passed from the MABL model back to the atmospheric 
model for use in the next time step. Fig. 2.6 shows the 
interrelationships between the atmospheric and oceanic 
models in the coupled model. 

The coupling process has the potential of significantly 
enhancing the OBL prediction qualities because of the 
improvement in the boundary conditions inherent in a coupled 
OBL-MABL system. Longer period reasonable forecasts should 
be obtained due to the feedback between the adjacent turbu- 
lent boundary layers and the associated thermodynamic and 
dynamic adjustments in each layer. Formation of stratus or 
fog in the MABL model is particularly important. Accurately 



32 



r 







Fig. 2.6 



Flow Chart for Ccupled OBL and MABL Model 
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E. 0* ESTIMATION FBOM TURBOIENT KINETIC ENERGY (TKE) 
DISSIPATION 

Direct measurement of u* on a ship is rather difficult. 
One popular technique which is used extensively at the Naval 
Postgraduate School is the dissipation method. An excellent 
review of this method as well as the bulk method is 
contained in Schacher, et al. (1981) . The dissipation 
method involves the use of the turbulent kinetic energy 
(TKE) equation which, assuming steady conditions and the 
divergence term to be small relative to the ether terms, 
becomes ; 

u * = {£€kz ]/[ <t) m (z/L) - (z/L) ]}0 (2.37) 

The denominator on the right hand side of equation 2.37 has 
been parameterized, where the quantity (J) m (z/L) - (z/L) is 

substituted by following Nyngaard and Cote (1971): 

(J) 6 = (1 + 0.5 Iz/Lj^) 3 / 2 for z/L < 0 (2.38) 

(J> € = (1 + 2.5 (z/I) ^) ^ for z/L > 0 (2.39) 

Determination of the dissipation rate, €, necessary for this 
technique requires that the inertial subrange follows a 
Kolmogorov spectrum, i.e., where: 

s u (k*) = a e ^k*-5/3 (2.40) 

S u (k*) in equation 2.40 is the spectral density of the hori- 
zontal windspeed, <* is a coefficient with a magnitude of 
0.52 (Champagne, et al- , 1977), and k* is the wavenumber. 

Equation 2.40 has been experimentally verified. On the R/V 
Acania a hot film sensor was used to obtain a spectrum of 
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high frequency velocity fluctuations from which the TKE 
dissipation u* was calculated. It is important to note here 
that poor velocity turbulence spectra are obtained during 
periods cf fog, rain cr drizzle. This leads to the dissipa- 
tion method being invalid. Errors may also occur if the 
relative wind is less than 3 m/s allowing condensation to 
form on the hot film. 

Applying Taylor’s hypothesis, i.e., k* = 2*f/o, equation 
2.40 may be rearranged to become: 

€ = 240 [ S u (f) ]0 (f/0) 0 (2.41) 

where f is frequency measured iD Hz. The inertial subrange 
spectral densities are integrated between two wavenumbers, k^ 
and k h respectively, in equation 2.40 following Khalsa and 
Bussinger (1977) : 

/ k h 

s u (k)dk = (3/2) a € 2/3(k/2/J - k-0) (2.42) 

Applying Taylor’s hypothesis, equation 2.42 now becomes: 

a Af 2 = 0. 230 (€?) 0 (f/0 - f"0) (2.43) 

Rearranging equation 2.43 and combining with equations 2.37, 
2.38 and 2.39, with a von Karman constant of 0.4, we obtain: 

u* = 2.81 a u (z/[ U<|) (z/L) ]} 0 {f/0 - f/0} 1 / 2 (2.44) 

Substituting f = 50Hz, f = 5Hz, the measurement height, z = 
20.0m, equation 2.44 reduces to its operational form as 
follows: 

u* = 3.69 a [ U 20m (t € (z/L) J"0 (2.45) 
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Because z/L is a function of the friction velocity, u*, 
according to equations 2.38 and 2.39, the solution to equa- 
tion 2.45 involves an iterative process {Large, 1979). 
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III. SYNOP TIC SITUATION 

A. MIXED LAYER DYNAMICS EXPERIMENT (MILDEX) 

The Mixed Layer Dynamics Experiment (MILDEX) was 
conducted during the period of 24 October to 10 November, 
1983 near 34°N and 126°fl. MILDEX was a multi-group and 
multi-platform experiment designed to: 

• provide magnitudes of air-sea energy exchange rates for 
use as boundary conditions in mixed layer modelling. 

• evaluate the drag coefficient from TKE dissipation meas- 
urements and its dependence on swell amplitude and 
direction. 

• to provide a time series, from radiosonde measurements, 
of MABL structure for model verification. 

• provide an intercomparison between radiative transfer 

algorithms developed at Scripps Institute of 

Oceanography and measured data. 

• evaluate the utility of SODAR as a ship-borne 

instrument. 

• evaluate radiative transfer models of the atmosphere 
with cloudiness as a unique parameter. 

Meteorology data were collected on the R/V Acania as well as 
two other platforms during MILDEX. R/V Acania data collec- 
tion consisted of measurements of windspeed, temperature, 
humidity, radiation and atmospheric pressure in the surface 
layer, and sea surface temperature (measured by thermistor) , 
and radiosonde derived profiles of atmospheric windsjeeds 
and directions, temperatures and humidity. Radiosonde 
launches were coordinated with satellite pass times. 
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Fig. 3. 1 depicts the locations on the bow of the R/V 
Acania of mean and turbulent measurement sensors. There 
were two levels of these measurement sensors: bow mast at 5 
meters and main mast at 20 meters. In addition, there were 
hourly observations of sea state and cloudiness. 
Conductivity, temperature and depth (CTD) casts were also 
made from the R/V Acania. Two ether vessels, the R/V Wecoma 
and R/P FLIP collected meteorological and oceanic data near 
the same site as R/V Acania. During the periods of interest 
in this thesis, the R/P FLIP was always within two kilo- 
meters of the R/V Acania. The meteorological data has been 
compiled in the MILDEX report (Geernaert, et al. , 1983) . 



MET. TOWERS 

MEAN AND 
TURBULENCE 
SENSORS 




Fig. 3.1 Schematic of the R/V Acania 
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Time series of wind speeds, directions, SST, air temper- 
ature, and surface pressure in Fig. 3.2 show seme of the 
variability experienced during MILDEX. Frontal passages are 
evident on 29 October, 31 October, and 06 November. For 
most of the experiment, the SST was greater than the air 
temperature (except cn 29 October and 09 November) , so the 
surface layer was unstable or near neutral- Only for 6-8 
hour periods on 29 October and 09 November was the surface 
layer stable. 

Three periods (cases) were considered for purposes of 
this study. These cases were chosen on the basis of strong 
wind forcing of the OBL and the availability of reliable 
radiosonde data for model initialization. Fig. 3.3 shows 
the movement of the R/V Acania for all three cases. The R/V 
Acania was located near 34. ON, 125. 4W during all periods. 
All times a re Pacific Dayli ght Time unles s otherwise speci - 
fied . MABL turbulent kinetic energy dissipation rate 
derived u*’s were obtained from velocity fluctuation meas- 
urements from the R/V Acania main mast (at 20 m) and a bow 
mast (at 5 m) instruments. MLB data were obtained from CTD 
casts from the R/V Acania for Case I and by the R/P FLIP for 
Cases II and III. 

Conditions associated with these cases was dominated by 
increasing clouds, swell and winds with intermittent to 
heavy rainfall. Heavy rain showers occurred on 29 October 
when the first of several fronts passed through the area. 
Since the presence of fog, rain or drizzle affects the 
validity of the hot film signals used in the dissipation 
technique, care was taken to exclude affected results during 
the periods. 

Ocean conditions were also variable. MILDEX experiment 
participants reported large scale advection and upwelling 
occurred during the course of the experiment. Significant 
internal wave activity was also evident and is being 
analyzed by other investiga tors . 
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Fig. 3.2 Time Series of MILD El Hind Speed and Direction, 
SST, Air Temperature, and Surface Pressure. Times are PDT. 
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Fig. 3.3 Movement of R/Y Acania: Case I (solid) , 

Case II (dash) , Case III (dot) . 

B. CASE I (0900, 31 OCTOBER - 0900, 01 NOVEMBER) 

During this period, the R/V Acania experienced weak ccld 
frontal passage at about 1700 on 31 October. Fig. 3.4 shows 
the National Meteorological Center (NMC) surface analysis 
for Case I. The synoptic situation was that a low center 



was located in the northeast Pacific ocean, and a high domi- 
nated the Pacific north of Hawaii. The entire period had 
nearly 100/? cloud cover and intermittent rain or drizzle. 
Several rain squalls were noted after the frontal passage. 
Clouds varied from altostratus prior to the frontal passage 
to stratus for the rest of the Case I period. No fog was 
reported, so a considerable amount of hot film data was 
available. Fig. 3.5 is a copy of GOES WEST visual imagery. 

As shown in Fig. 3.2, winds were initially from the SSW 
with a shift to the west after the frontal passage with 
speeds varying from 7 to 9 m/s ahead of the front and from 3 
to 5 m/s behind the front. There was a steady decrease in 
surface pressure from 1018 mb tc 1014 mb. Swells increased 
from 4 to 9 ft after the frontal passage and gradually 
decreased to 5 ft towards the end of the observational 
period. Air temperature dropped by 1.5°C after the frontal 
passage. Throughout the period, the SST was greater than 
the air temperature with the greatest difference at night as 
expected. Air-sea temperature differences varied from -.2° 
to -2.3° C. The MABL surface layer was near neutral to 
unstable during the period. 

C. CASE II (160 0, 03 NOVEMBER - 1600, 04 NOVEMBER) 

During this period, a stronger cold front than in Case I 
was approaching. Fig. 3.6 is the NMC surface analysis for 
Case II. A low center was positioned in the Gulf of Alaska, 
and a high pressure center had intensified west of Baja 
California and south of the Acania. No frontal passage was 
observed for this case. Winds were generally from the NNW 
shifting to the west in the middle of the observational 
period and then returned to the NNW at the end while wind 
speeds fluctuated between 3 and 7 m/s. The surface pressure 
was very steady changing only by 2 mb from 1021 to 1023 mb. 
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Swells were significant and from the NW and ranged from 7 to 
9 ft. Fluctuations in air temperature never exceeded 1.0°C. 
As in Case I, the SST was greater than the air temperature 
with the air-sea temperature difference varying from -.6°to 
-1.6°C. In this case, the MAEL surface layer was unstable 
throughout the entire period. 

Fig. 3.7 shows GOES WEST visual imagery for this case. 
Cloud cover was heavy in the late afternoon and again at 
mid-morning on the next day. Curing the night, the sky was 
cloud-free. Fog, cumulus and stratus clouds were evident in 
the afternoon. Stratus with some cirrus was present the 
next morning. No precipitation was recorded during the Case 
II period. 

D. CASE III (1600, 04 NOVEMBER - 1600, 05 NOVEMBER) 

The cold front approaching during Case II dissipated at 
the beginning of the period and the area came under the 
influence of a weak high pressure south of the R/V Acania. 
Fig. 3.8 is the NMC surface analysis for Case III. Another 
low pressure system is to the NW with an associated cold 
front. This second system did not appear to influence the 
weather or the seas in this case. Winds initially were from 
the NNW and gradually shifted to the SSW as the high moved 
through the area. Speeds generally decreased from about 7 
to 2 m/s, and air temperature fluctuated no more than .7°C. 
Surface pressure varied by 2 mb from 1022 to 1024 mb while 
sea swells were from the NW at 8 to 10 ft. As in the 
previous two cases, the SST was greater than the air temper- 
ature with the air-sea temperature difference varying from 

O O 

-1.1 to -2.1 C. Again, the MABI was unstable for the entire 
period. 

This last case had the least cloud cover due to the 
influence of the high pressure system. Fig. 3. 9 shows GOES 
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WEST visual imagery for Case III. At the beginning of the 
period, 50% cloud cover, consisting mainly of stratocumulus, 
was evident. This changed to cirrus and finally a clear 
night-time sky. On the morning of 5 Nov. there was 30% to 
50% cloud cover (mostly cirrus and some cumulus) with 
clearing in the afternoon. Nc precipitation was recorded 
for this period which was good for hot film performance. 
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Fig. 3.4 NMC Surface Analysis 1200 GMT OCTOBER 31,1983 
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Fig. 3.5 GOES kEST Visual 2215 GMT NOVEMBER 01, 1983 






47 



Fig. 3.6 NMC Surface Analysis 1200 GMT NOVEMBER 03, 1983 
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Fig. 3.7 GOES SEST Visual 2215 GMT NOVEMBER 03, 1983 
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Fig. 3.8 NMC Surface Analysis 1200 GMT NOVEMBER 04, 1983 
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GOES IEST Visual 2315 GMT NOVEMBER 04, 1983 












17. DATA AN I RESOLTS 

A. APPROACH TO FORCING 0* INTERPRETATION 

For each of the three cases, the u* associated with the 
model predictions was compared with the TKE dissipation rate 
(e) derived u* and the model MLD was compared with the 
observed MID. As described previously, the epsilon (€) 
estimates necessary fcr u* calculations by this method were 
obtained from hot film measurements on the bow (5 m) and 
main masts (20 m) of the R/V Acania (Fig. 3. 1) . 

The MABL model is configured for ten input winds, so the 
ten greatest and least winds were chosen from each 24-hcur 
period. These will be referred to in the discussion as 
"peak” winds. Linear interpolation was used for data points 
between these ten points. However, a forecaster at sea 
would not predict high and low winds, but would most likely 
forecast average winds over the period. Therefore, another 
model run utilizing 2.5 hour wind averaging was input over 
the same periods to see if the use of averaged winds would 
have any deleterious effect on the MLD computations. 

The purpose of these comparisons is to evaluate atmos- 
pheric forcing (u*) in relation to predicted and observed 
MLD 1 s . There are three u*'s used in this study: those 
estimated by bulk formulations using peak winds, by bulk 
formulations using averaged winds, and by dissipation esti- 
mates. In the following discussion, the first two sets will 
be referred to as model u* (peak or averaged) and the last 
set as dissipation u*. 



B. CASE I (09 00 , 31 OCTOBER - 0900, 01 NOVEMBER) 

This first case resulted in the poorest comparison of 
predicted and observed values of any of the three cases. 
Initialization data is summarized in Table I. Table II 
lists the 2.5 hour averaged winds that were used for subseq- 
uent runs. In Fig. 4.1, a comparison of the two wind inputs 
is shown. Fig. 4.2 shows a comparison of u**s. A diffi- 
culty in comparing model and instrument derived values for 
this case was instrument performance and equipment break- 
down. Both of these problems were encountered. 
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TABLE I 

Case One Initialization 



Atmospheric Val ues 

Pressure: 1019.0 mb Sea Temp (C) : 18.64 

Pot. T (C) : 17.9 Spec. Hum: 10.47 

Lift. Cond. Lvl: 351.74m Aavection: 0. 

Oce an Values 

Mixed Layer Depth: 34.5m Temp. Jump (C) : -.16 
Below Layer Gradient: -.0155 

Wind s 

Time (PDT) S pe ed (m/s) D irect ion 



0900 


8.3 


1 120 


9.5 


1345 


7.4 


1745 


9. 1 


2030 


3.4 


2130 


6.4 


0140 


3.5 


0345 


6.1 


0505 


2.4 


0720 


5.9 



181.0 

176.0 
205. 0 

199.0 

218.0 

215.0 
259. 
247. 
230. 
294. 



0 

0 

0 

0 



Rain, drizzle, and low relative wind speeds affected the 
accuracy of both the bow and main mast dissipation derived 
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Wind Speed (M/S) 













TABLE II 








Case One Averaged 


Hinds 




lime (PDT) 


Speed (m^s) 


Direction 




0900 


8. 3 


181.0 




1130 


8.9 


178.0 




1400 


8.2 


203. 0 




1630 


8.6 


196. 0 




1900 


4.6 


263.0 




2130 


5.4 


221.0 




0000 


5.7 


253. 0 




0230 


4.5 


253. 0 




0500 


4.2 


231.0 




0730 


4.9 


295.0 






Fig. 4.1 Comparison of Case I 
Initial and Averaged Hinds. 
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Fig. 4.2 Case I Friction Velocities 

u* values. Post-processing removal of the bad data resulted 
in significant data gaps. For the bow mast, dissipation u* 
values are missing from 1616 tc 2315 and from 0604 onward. 
Only 12 useable observed data points were available (as 
compared to 48 for the model) for the period. The field log 
indicates a lack of confidence in main mast spectra by the 
observers onboard due to spurious voltage readings and 
connector problems, so data from it were not included at all 
in this comparison. 

Fig. 4.2 shows the peak wind and averaged wind model u* 
and the dissipation u* values. As expected, the model u* 
values are a function of the differences in the input winds. 
Where the averaged winds are higher, the u*'s are higher. 
The dissipation u* values obtained from the bow mast have 
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the same trend but are lower than both model u*'s in the 
pre-frontal period (before 17 0 0) . They are close to the 
"peaked” wind model predictions for the rest of the time. 
All of the u* dissipation values are lower than the averaged 
wind model predictions. 

These results of model and dissipation u* values were 
not completely expected for the frontal regime. Previous 
comparisons of bulk and dissipation u*’s in frontal regions 
(Geernaert, et al. , 1985) indicated that the dissipation 
u*'s were higher than bulk values prior to a frontal passage 
and then the magnitudes reversed after the passage. Here 
the dissipation u*’s remained lower than the bulk value 
throughout. Generally, the study is revealing that small 
pre-frontal wind speeds (less than 6 m/s) produce greater u* 
deviations (turbulent vs bulk) . Additionally, if the swell 
and wind oppose each other, greater u* deviations were also 
evident. These were minimized if the wind and swell were at 
right angles to each other. The reasons for this are not 
entirely clear. However, it has been postulated (Geenaert, 
et ai., 1985) that the direction and magnitude of the swell 
relative to the wind (primarily the interaction of wind 
stress with short wave elements on long wave crests) affects 
the roughness length, z o , which in turn affects u* and C d . 

No attempt will be made with these results to examine 
how the swell interacted with the wind and how this could 
have modified the u* values. Perhaps this case did not have 
the same kind of swell-wind interaction as in previous 
studies. This is an area for further investigation. 

Fig. 4.3 compares the model and observed HLD’s. The 
model prediction was based on bulk u* values shown in Fig. 
4.2. In this case, model results (solid line) early in the 
period do not agree with the observed MLD (dashed line) . 
While both curves nearly coincide at the start, the model 
curve shallows much more quickly although it reaches the 
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Fig. 4.3 Case I Mixed Layer Depth 

proper depth. Then the model deepens rapidly and stays 
significantly deeper for the rest of the case. An MID error 
of up to 10 meters can be tactically significant for naval 
forces in predicting sonar performance. 

The early shallowing may not constitute a model 
prediction error. A CTD cast would give the thermal MLD 
while the model predicts the turbulent boundary layer depth. 
It takes a finite amount of tiae for the temperature field 
to adjust to the turbulent field. The expression for the 
lag is: 

t ~ [ (aT) u** ]/Q 0 2 (4.1) 

where : 

T = temperature 
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t = lag time 

Q o = surface heating 

Thus based on equation 4. 1, the greater the u*, the greater 
the lag and the greater the heating (Q o ) , the less the lag. 

Using 2.5 hour averaged winds showed no improvement in 
MLD prediction (Fig. 4.3 dotted line) and- in fact gave a 
poorer prediction. While the two wind inputs (Fig. 4.1) are 
quite different, the output MIL’S are quite close together 
except at 1300. This is a pre-frontal period of strong 
winds and, from 1200 to 1700, the averaged wind speeds are 
higher than the initial winds. Since MLD is a function of 
the wind speed cubed, the stronger winds would give deeper 
mixing. This is again true for the period from 2300 to the 
end of the run where the averaged winds are generally higher 
than the initial winds. Overall, one can conclude that for 
this case both the "peaked" and averaged winds give poor MLD 
predictions and these could have adverse effects for naval 
operations. However, if the dissipation u* values for 0900 
to 1700 (which are lower than the bulk values) had been used 
by the model, the model MLD would have been closer to the 
observed. 

C. CASE II (1600, 03 NOVEMBER - 1600, 04 NOVEMBER) 

The coupled model MLD prediction agreed well with the 
observed MID changes in this case. The model predicted 
unstable ABL conditions throughout the period which agreed 
well with observations. Table III contains the initializa- 
tion data for this case and the 2.5 hour averaged winds are 
listed in Table IV. Fig. 4.4 compares the two input winds 
for this case. 

This was a period with unreliable dissipation u* data. 
Fig. 4.5 is a plot of the various friction velocities. 
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TABLE III 








Case Two Initialization 






Atmospheric Values 




Pressure: 1022.4 mb Sea 


Temp (C): 18.78 
. Hum: 9.32 




Pot. T (C) : 
lift. cond. 


17. 1 Spec 




Lvl: 476.71m Advection: 0. 






Ocean Values 






Mixed Layer 


Depth: 12.0m Temp 


. J ump (C) : - . 2 




Below Layer 


Gradient: -.0094 






Winds 






Time (PDT) Speed (m/s) 


Direction 




1600 


4.0 


040.0 




1745 


5.6 


310. 0 




2120 


8.6 


036. 0 




2230 


4.7 


319.0 




2300 


6.4 


355.0 




0135 


2. 4 


270. 0 




0400 


7.0 


281.0 




0705 


3.0 


350. 0 




1005 


7.2 


008. 0 




1420 


4.6 


336. 0 





While fog or rain were not a factor in dissipation measure- 
ments, equipment failure was a problem. This led to a 
reduction of dissipation u* data. The bow mast data were 
not available from 1700 to 2300, 2330 to 0217, and 0900 to 
1600 when a major equipment repair was accomplished. Bow 
dissipation u* data from 12 thirty minute periods were 
available. Main mast dissipation measurements were not 
available from 1600 tc 0752 and from 0822 to 1340. It was 
also absent after 1410 and only data from 4 thirty minute 
periods were available. 

At the beginning of the period, the bow mast dissipation 
u**s are close to the bulk derived (model) values. However, 
during the middle part of the period, the bow dissipation 
values are consistently lower than the bulk (model) although 
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lime (PDT ) 


TABLE 17 

Case Two Averaged 
§.£S§d (m/s) 


Hinds 

Directi on 


' 


1600 


3.9 


040.0 




1830 


4.5 


345. 0 




2100 


5.7 


349.0 




2330 


4.3 


308. 0 




0200 


3.6 


274.0 




0430 


6. 1 


282.0 




0700 


4.6 


307. 0 




0930 


6.6 


357. 0 




1200 


6.8 


343. 0 




1430 


5. 4 


341. 0 






Fig. 4.4 Comparison of Case II 
Initial and Averaged Hinds. 
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Fig. 4.5 Case II friction Velocity 

they tend to have the same time variations. The two bcw u* 
dissipation values near 0800 are close to the model values. 

It is apparent from Fig. 4.6 that the model predicted 
the MID well enough for Naval tactical applications. The 
observed MLD (dashed line) shows some differences from 0200 
to 1000 which are probably due to internal wave activity, 
advection or as a result of the ship moving to a different 
area. Observed MLD data was not available for 1100 to 1300. 
This contributes to the abrupt MLD rise at 1400. The model 
(solid line) predicted the MLD guite well especially for the 
deepening part of the period. Better agreement between 
observed and model shallowing would probably be evident if 
MLD data would have been available for 1100 to 1300. 
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Fig. 4.6 Case II Eixed Layer Depth 



Comparison of the peaked winds with the 2.5 hour aver- 
aged wind MLD's yields some interesting results. During the 
deepening phase, the averaged winds consistently give a 
shallower MLD (dotted line) than the initial winds. At this 
time the "peaked" winds are consistently stronger (Fig. 4.4) 
than the averaged winds (except at 0100) and would produce 
more vigorous mixing. However, during the shallowing 
portion, the averaged wind prediction is significantly 
deeper than the peaked wind prediction. This time the aver- 
aged winds are predominantly stronger than the peaked winds. 
Again the predicted MLD results are expected since the MLD 
is proportional to u* 3 . If the dissipation u* values for 
0200 - 0600 would have been used by the model, the predicted 
MLD would have been shallower fcr that period and closer to 
the observed. 



D. CASE III (1600, 04 NOVEMBEB - 1600, 05 NOVEMBEB) 



As in the previous case, the coupled model MID 
predictions agreed well with the observed conditions. ABL 
conditions were unstable throughout the case. Table V lists 
the initialization data utilized for this model run, and 
Table VI lists the 2.5 hour averaged winds. Fig. 4.7 
compares the two input winds. 





TABLE V 




Case Three Initialization 




Atmospheric Values 


Pressure; 1022.3 


mb Sea Temp (C) : 18.70 


Pot. T (C) : 17.7 

Lift. cond. Lvl: 


Spec. 


Ham: 9.48 


508.0m Aavection: 0. 




Ocean Values 




Mixed Layer Depth; 11.5m Temp. 

Below layer Gradient: -.0053 


J ump (C) : - . 1 




Winds 




Time (PDT) 


Speed (m/s) 


Direction 


1600 


5.2 


340.0 


1755 


6.5 


338.0 


2100 


5.6 


351.0 


2300 


4.5 


352.0 


0046 


3.7 


338.0 


0155 


5.2 


325.0 


0445 


2.6 


020.0 


0715 


2. 1 


296.0 


0915 


4.8 


333.0 


1045 


2.4 


271.0 

... 


This case had 


no fog related 


problems affecting the 



dissipation data but for some periods, the low relative wind 



speeds degraded dissipation reliability. Several large 
dissipation data gaps resulted. For the bow mast, data gaps 
exist for 1600 to 0018, 0100 to 0815, 0845 to 1200, and from 
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TABLE VI 






Case 


Three Averaged 


Binds 




TimejPDT) 


S£eedjm/s) 


Direction 
— 3IT0TU 




1830 


6*0 


340. 0 




2100 


5.6 


351.0 




2330 


4.3 


346. 0 




0200 


5.0 


3 54. 0 




0430 


2.7 


004.0 




0700 


2.3 


311.0 




0930 


3.2 


272. 0 




1200 


2.5 


264.0 




1430 


2.3 


235. 0 






Fig. 4.7 Comparison of Case III 
Initial and Averaged Binds. 
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Fig. 4.8 Case III Friction Velocities 



1331 onward. Only seven reliable dissipation values were 
calculated for the bow mast. Data was not available for the 
main mast from 1600 to 0018 / 0100 to 0315, and from 1410 to 
the end of the period. This resulted in four thirty-minute 
dissipation averages. Fig. 4.8 compares the different fric- 
tion velocities. 

For the entire period, the peak and averaged wind model 
u*'s are very close due to the small wind differences 
between the two (Fig. 4.7). Both the main mast and bow mast 
dissipation u* ' s agree reasonably with u*'s from the peaked 
wind model run. This is particularly true in the second half 
of the period. Some variation is expected since the dissi- 
pation u* comes from instantaneous measurements, while the 
model uses interpolated values between the input winds. 
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Fig. 4.9 Case III Hixed Layer Depth 

The model also predicted the HLD's well (Fig. 4.9) 
althouyh there is greater deviation from the observed 
(dashed line) than in Case II. As before, the actual SLD 
shows some internal wave or advective features (particularly 
around 0000 and 0400). Good observed MLD data was available 
for the entire period. 

The initial model run (solid line) agrees well in the 
initial stages but deviates during the above mentioned 
internal wave or advective events. The sharp "dip" at about 
0930, not present in the model output, is an artifact intro- 
duced by use of a spline interpolation in the plotting 
routine. The model also shallows to the proper depth but it 
occurs about one hour early. 



6 5 



Comparing the averaged winds ML D with the peak winds MID 
shows that both are very similar. The plot of input winds 
(Fig. 4.7) shows that both the initial and averaged winds 
were very close throughout the period except between 0800 
and 1000, so both model runs should produce similar MID’s. 
After 1100, the averaged wind MLD (Fig. 4.9 dotted line) 
shows a double maxima during shallowing. It is unknown at 
this time why the model predicted such maxima. They have 
been slightly exaggerated by use of the spline interpolation 
but are present in the model output. As with the initial 
wind MID, an anomolous "dip 1 * is observed at 1030 also due to 
the spline interpolation. Both the initial and averaged 
wind MID's are within operatioral tolerances for Navy use 
but the forecaster should be well versed in the model prop- 
erties and limitations. 
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7. CONC1QSIONS AND R ECO MMENDAT IONS 
A. CCNC10S1O1S 

The coupled model has been subject to various past anal- 
yses and performed well. This study has examined the 
performance of the model in regimes that stretch its capa- 
bilities. Several important conclusions can be drawn: 

(1) The model overpredicts clouds which affects short and 
long wave radiative effects on boundary layer 
dynamics. 

(2) Good comparisons of bulk and dissipation u* for 
unstable conditions are made. 

(3) In near neutral or slightly unstable regimes, u* 
predictions vary systematically probably due to swell 
and wind stress interaction. 

(4) Better predictions are achieved if peak winds instead 
of averaged winds are used. However, a forecaster 
would only have average wind prediction estimates 
available at the time of operationally initializing 
the model. 

(5) Given reasonable fluxes cf short and long wave radia- 
tion, the model predicts routinely observed varia- 
tions in ocean mixed layer depth. 

(6) It is a potentially powerful tool and an aid in 
single-station forecasting to a Naval Oceanographer 
who is aware cf the model limitations. 

The model performed very well within its limitations. 
However, further work needs to be done before it can be 
implemented in an operational role. 
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B. RECOMMENDATIONS 



In order to facilitate operational usage, several areas 
should be explored. Some of the aforementioned problems 
could be alleviated by attempting the following 
recommendations: 

(1) Continue research and development of the coupled 

model with direction toward operational 

implementation. 

(2) When clouds form, the model assumes them to extent 
from horizon tc horizon. This is unrealistic and a 
statistical probability of percentage cloud cover 
should be added. Partial cloud cover would also 
alter short and long wave effects on boundary layer 
dynamics. 

(3) Incorporate swell and wind interaction effects into 
u* calculations based on bulk formulae. 

(4) Continue to analyze MILEEX data, particularly for 
periods where there are more complete dissipation 
u*'s. 

(5) Plan for model coupling to other predictive models 
such as IREPS and ICAPS. 

(6) Involve students in at-sea data gathering experi- 
ments. This would enhance understanding of thesis 
research and minimize "false starts" due to 
attempting cases with faulty or limited data. 

The model has tremendous potential. Performing further 
analysis and tweaking the model will result in a powerful 
tool for the fleet; one that can be used in times of crisis 
when extensive satellite networks may not be available for 
use by forecasters. 
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